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ABSTRACT 

Optically thin and geometrically thick accretion flows are known to be re- 
sponsible for the observed radio/X-ray luminosity relation of the X-ray bright 
galactic nuclei. It has also been suggested that supermassive black hole masses 
can be estimated from measurements of the core radio luminosity and the X-ray 
luminosity by using the advection-dominated accretion flow (ADAF) model. In 
this study we increase the number of data available by compiling the radio/X- 
ray fluxes and the mass in published literatures, and compare the observed ratio 
of the luminosities with predictions from various models of optically thin accre- 
tion flows. Semi-analytically derived relations of the luminosities are presented 
in cases of the standard ADAF model and modified ADAF models, in which 
a truncation of inner parts of the flows and winds causing a reduction of the 
infalling matter are included. We show that the observed relation can be used 
indeed to estimate the supermassive black hole mass, provided that properties of 
such accretion flows are known. Having investigated sensitivities of the method 
on modifications of the 'standard' ADAF model, we find that a general trend 
of model predictions from the 'standard' ADAF, the truncated ADAF and the 
'windy' ADAF are somewhat indistinguishable. We also find, however, that the 
extreme case of the windy model is inconsistent with currently available obser- 
vational data, unless microphysics parameters are to be substantially changed. 
High resolution radio observations, however, are required to avoid the contamina- 
tion of non-disk components, such as, a jet component, which, otherwise, results 
in the over-estimated SMBH mass. 

Subject headings: accretion disks - black hole physics - galaxies:nuclei - 
radiation mechanisms:nonthermal - radio continuum:galaxies - X-rays:galaxies 
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1. INTRODUCTION 

Supermassive black holes (SMBHs) have been considered as the most likely power 
sources of the activity in quasars and Active Galactic Nuclei (Lynden-Bell 1969; Rees 
1984). SMBHs at the centers of all galaxies are now recognized as ubiquitous, whose mass 
^fsMBH is proportional to the spheroidal bulge mass of the host galaxy or the galactic bulge 
luminosity (Kormendy & Richstone 1995; Magorrian et al. 1998; Richstone et al. 1998) 
and is strongly correlated with the velocity dispersion of the host galaxy (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000a; Ferrarese 2002). Recently, evidence for the existence 
of a SMBH in the center of our Galaxy has been added (Eckart & Genzel 1997; Genzel et 
al. 1997; Ghez et al. 1998). Though searches for SMBHs are primarily based on spatially 
resolved kinematics, another possible way to infer the presence of SMBHs and to shed 
light on the physical condition is to examine spectral energy distribution over a wide range 
from the radio to the hard X-ray frequencies. This emission spectrum is produced by an 
accreting matter, as the surrounding gas accretes onto the central SMBH (e.g., Frank et al. 
1992; Ho 1999). 

Several attempts have been made in utilizing this fact to estimate SMBH masses (Yi & 
Boughn 1998, 1999; Franceschini et al. 1998; Lacy et al. 2001; Ho 2002; Alonso-Herrero et 
al. 2002; Cao & Jiang 2002). It has been also suggested that advection-dominated accretion 
flows (ADAFs) are responsible for the observed radio and X-ray luminosities of some of the 
X-ray bright galactic nuclei (Fabian & Rees 1995; Di Matteo & Fabian 1997; Yi & Boughn 
1998, 1999). Moreover, Yi & Boughn (1999) discussed that the observed radio/X-ray 
luminosity relation in terms of the 'standard' ADAF model and demonstrated that it could 
be used as a relatively effective and consequently suitable tool to estimate unknown black 
hole masses. In this study we further explore this possibility both by increasing the number 
of data available and by deriving similar relations for the models with additional features, 
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such as, truncations, winds, and investigate sensitivities of the method Yi & Boughn 
suggested on modifications of the 'standard' ADAF model. 

This paper begins with model descriptions for the standard ADAF model and modified 
versions of the model in § 2. We derive analytical expressions to describe the radio/X-ray 
luminosities for various models. We present the calculated ratios of the luminosities and 
compare them with the observational data in § 3. We discuss and conclude in § 4. 

2. RADIO/X-RAY EMISSION FROM ACCRETING SMBHs 

When a mass accretion rate is below a critical rate, the radiative cooling becomes 
slower than the viscous heating. As a result of this, the dissipated accretion energy is 
inefficiently radiated away and advected inward to the central object with the accreted 
matter, which is a key idea of ADAF models (see Rees et al. 1982; Narayan & Yi 1994, 
1995a,b). Detailed numerical calculations based on this idea have been performed, and the 
resulting spectra have been successfully applied to a number of extra galactic systems (e.g., 
Narayan et al. 1995; Lasota et al. 1996; Mahadevan et al. 1996; Manmoto et al. 1997; 
Narayan et al. 1998). 

In this study we apply analytical scaling laws for self-similar solutions to derive the 
luminosity relation for ADAF- variations (Narayan & Yi 1994; Mahadevan 1997). We adopt 
the following dimensionless variables : mass of the SMBH m = M/M Q , where M Q is the 
solar mass; radius from the SMBH r = R/R g , where R g = 2GM/c 2 = 2.95 x 10 5 m cm; 
and mass accretion rate m = M/M Edd , where M Edd = L Edd /r] cS c 2 = 1.39 x 10 18 m g s _1 
when 7] cS = 0.1. We assume that the flows are spherically symmetric as the ADAF 
is quasi-spherical (Narayan & Yi 1995a). However, the anisotropic extinction of host 
galaxies is not taken into account explicitly. This fact should be considered carefully if 
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one wishes to elaborate this method. The advection fraction / is determined such that 
the ion and electron energy balance equations are met. The quantities of interest are 
the volume-integrated quantities which are obtained by integrating the heating rate and 
cooling rates throughout the volume of the flows. For spherically symmetric flows the 
volume-integrated quantities are defined by 

Q x = AnR 2 q x dR 

r 2 q x dr, (1) 

min 

where X denotes any quantity of interest (cf. Mahadevan 1997). As canonical values in 
a model for the ADAFs parameters are taken to be r min = 3, r max = 10 3 , a = 0.3, and 
(5 = 0.5, where r min and r max represent the inner and the outer boundaries of the disk, 
respectively, a is the viscosity parameter, and j3 is the ratio of the gas pressure to the total 
pressure. The results of this study are insensitive to these values unless these parameters 
are very small. The lower limit of the integration corresponds to the last stable orbit for 
the Schwarzschild black hole, and the upper limit is not very crucial since most of radiation 
and heating occur around the central black hole. We assume that the electron temperature 
is constant for r < 10 3 , as suggested by Narayan & Yi (1995b). 



2.1. Emission from 'Standard' ADAFs 

The total heating by the viscous dissipation for ions Q + is given by 

Q+ = 9.39 x l0^^l C3 mmr-] n , (2) 

where c 3 is a constant defined as in Mahadevan (1997). We neglect the possibility of 
the viscous heating of the electrons, since we expect that the fraction of viscous energy 
transferred to the electrons is in the mass ratio of the electron to the ion, ~ 1/1800. The 
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electrons are assumed to gain energy only by Coulomb interactions with the ions. The total 
heating of the electrons due to the Coulomb interaction is given by Stepney & Guilbert 
(1983) 

Q ie = 1.2 x ltf & g{d e )a- 2 cfczf3mm 2 r^ (3) 

where g(6 e ) = (2 + 26 e + 6' 1 ) exp(-9- 1 )/K 2 (e- 1 ), 6 e = kT e /m e c 2 , K 2 being a modified 
Bessel function, k, m e , and c being the Boltzmann constant, the electron mass, and the 
speed of light, respectively, C\ is a constant given in Narayan & Yi (1994). 

The synchrotron photons are self-absorbed and give a blackbody radiation upto a 
critical frequency v c . For a given T e , the synchrotron spectrum L sync is given by 

L s J nc = s 3 ( SlS2 ) 8 / 5 m 6 / 5 m 4 / 5 T e 21 / 5 z/ 2 / 5 , (4) 

where 8l = 1.42 x lO 9 ^ 1 / 2 ^ - /3) 1/2 q 1/2 cJ /2 , s 2 = 1.19 x l(T 13 :r M , and s 3 = 1.05 x 1(T 24 , 
xm = 2z//3z/{,# 2 , i/ft = eB /2irm e c, e and S being the electron charge, the magnetic 
field strength, respectively. The radio luminosity L R at v is defined by uL s u ync . The 
highest radio frequency arises from the innermost radius of the accretion flows, r min ~ 3; 
u p = s\S2vrr x l 2 m x l 2 T'lr r ^^. At this peak frequency the peak radio luminosity is given by 

L R ee v v U^ = s 3 ( Sl s 2 )^/ 2 m^ 2 Tjr m ]i 4 . (5) 

We obtain the total power due to synchrotron radiation from 

Psync ^ P ' L s r c du=(5/7)u p L s y p nc . (6) 
Jo 

We set z/ min = 0, since z/ min is much smaller than v v . 

The bremsstrahlung emission per unit volume is given by Svensson (1982). 
Bremsstrahlung emission is due to both electron-electron and electron-ion interactions. The 
total bremsstrahlung power is given by 

P brem = 4.74 x 10 u a~ 2 cf\n(r max /r mm )F(e e )mm 2 , (7) 



-7- 



and the spectrum due to bremsstrahlung emission is given by 

L brem = 22Q x 10 24 a -2 c -2 ln(w / r min ) F (6 e )mm 2 T; 1 exp (-hu/ kT e ) , (8) 

where 

/Of) \ 1 / 2 

F(9 e ) = 4(-^J (l + 1.781^ 34 ) + 1.73^/ 2 (l + l.l^ + e e 2 -1.25^/ 2 ), e < 1, 

= ( ^J[ln(l. 123# e + 0.48) + 1.5] + 2.30 e (lnl.l230 e + 1.28). e > 1. (9) 

We neglect the Comptonization of bremsstrahlung emission and consider the 
Comptonization of synchrotron emission alone. A contribution by Compton up-scattered 
synchrotron photons to the hard X-ray luminosity becomes important as m increases, while 
bremsstrahlung emission dominates the hard X-ray luminosity when m is substantially 
low. We approximate the Comptonized spectrum by assuming that all the synchrotron 
photons to be Comptonized have an initial frequency of v v . The maximum final frequency 
of the Comptonized photon is Vf = 3kT e / h, which corresponds to the average energy of the 
photon for saturated Comptonization in the Wien regime. On average, we assume that all 
the photons would see one half the total optical depth (Mahadevan 1997), which is written 
as r es = 6.2a:~ 1 c]~ 1 mr~ 1 / 2 . The spectrum of the emerging photons at frequency v has the 
power-law shape 

L c v omp ~ L Ui {-)~ a % (10) 

where L Vi is the initial luminosity at the initial frequency Vi of the photons that are 
Comptonized, a c = — lnr es /lnA, A being the mean amplification factor in a single 
scattering. The total Compton power is therefore given by 



l>3kT e /h 

Pcomp / Li v ^ dv 

J u„ 



1 - a c \ 



6.2 x 10 10 T e 



l—a c 

11) 
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For given m, rh, a, (3, the total heating of the electrons should be balanced to the sum 
of the individual cooling; Q ie = P sync + Pbrem + P C om P - The electron temperature is obtained 
by adjusting until this equality is satisfied. 



2.2. Emission from Truncated ADAFs and 'Windy' ADAFs 

ADAFs have the positive Bernoulli constant, which is a kind of total energy of the 
fluid, and, therefore, are susceptible to outflows (e.g., Narayan & Yi 1994). It might be 
the case that jets are present near the inner region of ADAFs. Radio emission due to 
ADAFs at the high frequency will be suppressed if the inner part of the accretion flow is 
truncated by low frequency jets. In order to simulate this situation, we truncate ADAFs at 
a certain inner radius. We suppose that the truncation occurs at ~ 25 r g . An observation 
of NGC 4486 (M87) indicates that the jet is formed in a smaller radius than r R3 30 r g 
(Junor et al. 1999). This is also about where is the maximum radius that 15 GHz radio 
emission is generated for m = 10 7 ,m = 10~ 3 ,T e = IQ 9 K. Note that r can be estimated 
by r = (l/60)~ 4 / 5 x (i//15 GHz)" 4 / 5 (m/10 7 )- 2 / 5 (m/10- 3 ) 2 / 5 (T e /10 9 ) 8 / 5 . We consequently 
set r min in expressions in the previous subsection to 25 r g instead of 3 r g in corresponding 
expressions in the truncation case in order to obtain the radio luminosity and the X-ray 
luminosity. Since the frequency at which the radio luminosity is the highest arises from the 
innermost radius of the accretion disk, v p oc r min , the truncation radius affects the peak 
radio luminosity at this frequency according to Eq. (5). 

In addition to a jet or outflows, since the gas in ADAF solutions is generically 
unbound winds may also carry away infalling matter (see, e.g., Blandford & Begelman 
1999; Di Matteo et al. 1999; Quataert & Narayan 1999). We allow fn to vary with r as 
m = m out (r / r max ) p , where m out is the mass accretion rate at r max , p is the parameter 
which controls the effects of winds (see, e.g., Blandford & Begelman 1999). The mass loss 
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causes a significant effect on a model spectrum for a large value of p. Bremsstrahlung and 
synchrotron emissions decrease with increasing p, due both to the lower density and to the 
lower T e in inner parts. Compton emission decreases with increasing p even more strongly 
than the other two emissions since it depends both on v p L„ and a c . As m is replaced by 
riiout(r I r max ) p , expressions for the 'Standard ADAFs' are to be modified as described below. 
The total heating for ions Q + reads 

Q+ = 9.39xl0 38 ^c 3 mm Ott ,r-L^(C-i-r^), (12) 

where p < 1. Similarly, Q te is given by 

Q- = 1.2 x 10 38 ^ e )a-Vc 3 /?m<^ (13) 

where p ^ 1/2. We assume that xm and T e are constants with r and m as an approximation, 
as suggested by Narayan & Yi (1995b). Given x M , the cutoff frequency at each radius is 
given as v c = SiS2m~ 1 / 2 m^ t r m p JxT 2 r~^ 2p ~ 5 ^ A . The synchrotron spectrum is given by 

8 4p-6 -J— _*E_ 21 - 2 P ip-2 

L s u ync = s 3 (s 1 s 2 )^m^m 5 - t 2p r%a X 5 T e ^ 2p vW* . (14) 

The peak frequency and the radio luminosity at the peak frequency are correspondingly 
reduced so that we have 

VvKT = ^M'm^rhl^Tjr 6 ^ . (15) 
The total synchrotron power is also rewritten by 

Psync * ^Urdv= (^)^r- ( 16 ) 

For bremsstrahlung emission the total power and the spectrum are respectively given by 

P brem = 4.74 x l0 34 a- 2 cfmml ut r- 2 J> x F(6 e ) - r£j] , (17) 

L b J- = 2.29 x 10 24 a- 2 cr 2 mmLr-^F(^)T- 1 exp(-/u//fcT.) [^(r*. - rjj] . (18) 
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Note that as the exponent p approaches to 0, [^(r^ ax —r^ in )] becomes \n(r max /r min ). For the 

-v (J 2 p~ 1 )/ 2 _ 

max V max 



Comptonization of the synchrotron photons we take r es as r es = 6.2a c x rh out r, 
r mAn 1 ^ 2 )- The Compton spectrum and the total Compton power are respectively given by 



jComp 



comp 



V 



1 - a c \ 



6.2 x 10 10 T f 



1 — Ctc 



- 1 



(19) 
(20) 



3. BLACK HOLE MASSES AND RADIO/X-RAY LUMINOSITIES 

Figure 1 shows a plot of the ratio of the 15 GHz radio luminosity to the 2-10 keV X-ray 
luminosity versus the X-ray luminosity in logarithmic scales as the luminosities are derived 
in the last section. The solid lines are standard ADAF model predictions for SMBH masses 
of 10 6 — 1O 9 M , the dotted lines are ADAF model predictions with truncations at 25 r g , the 
short dashed lines and the long dashed lines are ADAF model predictions with winds, where 
p = 0.4 and p = 0.99, respectively. Values of p considered correspond to the moderate wind 
case, and the strong wind case. The case of p — 1 may represent the ADAFs where the 
convection is present in ADAFs in that the density scales —1/2 with the radius instead of a 
power of —3/2 (Narayan, Igumenshchev, & Abramowicz 2000; Quataert & Gruzinov 2000). 
For each curve, a constant fa for 'non-windy' models or m out for 'windy' models varies from 
10~ 4 (upper left corner) to 10~ L6 (lower right corner). A general trend of model predictions 
from the 'standard' ADAF, the truncated ADAF, and the 'moderately' windy ADAF is 
somewhat similar and overlapped. As long as the canonical model parameters remain 
unchanged, however, the 'strongly' windy ADAF model with large p, which represents the 
convection, results in quite different trends. Besides the trend, none of data points fall on 
the area defined by the convective ADAF model predictions. 

The open circles denote the observed core radio luminosity at 15 GHz, and the filled 
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circles the spatially resolved core radio luminosity which are converted from ones at 5 GHz 
to at 15 GHz using the z/ 7 / 5 power law (cf. Yi & Boughn 1999). We need such a conversion 
in order for a direct comparison of these luminosities at different frequency. Sources 
generally fall into the region where theoretical curves predict. For instance, the mass of the 
central black hole in NGC 3377 is known as about 8 x 1O 7 M (Kormendy & Richstone 
1995). According to our theoretical predictions from the models except the strong wind 
ADAF model, the mass of the central black hole in NGC 3377 falls in between 10 7 to 
10 8 M Q . Apparently, even though the inherent radio/X-ray luminosity relation provides an 
estimate of the central SMBH masses currently available, radio/X-ray flux relations are yet 
insufficiently accurate to give unambiguous estimates of the SMBH mass. Or, the observed 
radio/X-ray fluxes should be analyzed to extract 'disk components' more carefully. 

As for the sources which do not agree with the predictions very well, it is thought 
either because they have non-ADAF sources in their centers or because the observed radio 
and/or X-ray flux is due to a non-ADAF origin. The predicted mass of NGC 224 (M31) 
is smaller than the estimated mass at least by two orders of magnitude. NGC 224 has a 
very low core radio luminosity, which is even lower than predicted with the mass accretion 
rate that the X-ray luminosity implies on the basis of an ADAF model. NGC 224 seems to 
harbor a double nucleus (Lauer et al. 1993; Bacon et al. 1994). It is unclear whether ADAF 
is viable under such circumstance. Predicted masses of NGC 1068, NGC 3031 (M81), NGC 
3079, NGC 4151 are ~ 10 8 to ~ a few 10 9 M Q . On the other hand, reported mass estimates 
of SMBHs of these galaxies are less than ~ 10 7 M Q . Masses of these sources are thought 
to be over-estimated due to high radio/X-ray luminosities resulting from strong jet-related 
activities. Moreover, uncertainties in X-ray flux measurements and intrinsic variations could 
cause errors in our predictions. NGC 1068 is famous in that it has a strong compact radio 
source, which results in a much higher radio/X-ray luminosity ratio than other Seyferts 
(Gallimore et al. 1996). The radio core of NGC 3031 (M81) is also highly variable on many 
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timescale (Ho 1999). Moreover, NGC 3031 is substantially brighter in X-rays relative to 
the UV than in luminous AGNs (Ho 1999), and its core luminosity in the 2-10 keV band 
varies by a factor of ~ 2 in a timescale of years (Ishisaki et al. 1996). In NGC 3079 the 
high radio luminosity may be due to an unresolved small-scale jet. The radio morphology 
of NGC 4151 suggests that the emitting material has been energized from the nucleus and 
that the jet interacts with the ambient medium (Carral et al. 1990). 

In Figure 2, the predicted core radio luminosity is shown as a function of the SMBH 
mass in a logarithmic scale, with plots of SMBH candidates whose mass estimates are 
publicly known. From top to bottom, for each pair of lines, the mass accretion rate 
corresponds to 10~ 2 and 10~ 4 . The line types and symbols are same as those in Figure 
1. Although the radio flux data are in high angular resolution there appears a somewhat 
wide range of radio luminosities, which is likely the result of small, and thus invisible, 
radio jets or other components of various frequencies rather than the pure ADAF emission. 
In particular, the radio luminosity of objects such as NGC 3079, NGC 3031, and NGC 
1068 seems too high to be explained by the pure ADAF emission. In Figure 3, the X-ray 
luminosity for the 2-10 keV band is shown as a function of the SMBH mass in a similar way 
as in Figure 2. Quoted X-ray fluxes of NGC 4486 and NGC 3377 should be considered as 
upper limits (Reynolds et al. 1999; Pellegrini 1999). Note that the observed X-ray flux of 
NGC 3377 in the energy band of 0.2 - 4 keV is converted to the value in the energy band of 
2-10 keV, assuming a power-law with the canonical spectral index for Seyfert 1 galaxies, 
i.e., T = 1.7. High resolution X-ray images have revealed discrete X-ray sources for NGC 
224 (Trinchieri et al. 1999), NGC 3031 (Ishisaki et al. 1996), NGC 4736 (Roberts et al. 
1999). For these sources we adopt core or bulge X-ray fluxes. 
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4. DISCUSSION 

One of important predictions of the ADAF model is the radio/X-ray luminosity 
relation which can be used to estimate central SMBH masses. We compute the theoretical 
predictions of the radio/X-ray luminosity relation for the ADAF models with various 
modifications and compare the predicted relation with that from the available SMBH 
candidates. Several nearby extra galaxies have consistent ratios of radio to X-ray 
luminosities with the ADAF predictions for an estimated mass of the central SMBH. If the 
type of the accretion flow is provided by observations with a high angular resolution, the 
inherent radio/X-ray luminosity relation provides a direct estimate of the central SMBH 
mass. We show the logarithm of the ratio of radio to X-ray luminosities against the 
logarithm of the X-ray luminosity for several SMBH masses in Figure 1. The observations 
for several extra galactic objects are consistent with the predictions of ADAF models, 
as can be checked in Figures 2 and 3. It is, however, interesting to note that the model 
representing the ADAFs with the strong wind can be ruled out unless the microphysics in 
the standard ADAF model has to be modified significantly. 

In this study, we implicitly assume that all the SMBHs are non-rotating. When the 
black hole rotates, the last stable circular orbit is no longer 3 r g . It becomes smaller than 
this value and eventually the event horizon when the black hole spin parameter a ~ 1. If 
this is the case the boundary condition should change as well as the dynamical properties 
of the flows in inner parts (Gammie & Popham 1998; Popham & Gammie 1998). Besides 
dynamics of the flows, a rotating black hole may power a jet by the Blandford-Znajek 
process (Blandford & Znajek 1977). Therefore, even if the black hole rotates slowly, the 
rotation should be properly taken into account in order to make this method robust. 
Another crucial modification could be made in terms of microphysics we simplified. We 
have assumed that the viscously generated energy mainly heats the ions, heat exchanges 
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between ions and electrons occur only due to the Coulomb interaction, other physical 
processes, such as, MHD turbulence, pair production, are ignored. Many works have 
been done to quantify the effects of such processes (e.g., Ozel et al. 2000). Once clear 
understanding is incorporated in the method we present here, its implication can be more 
robust. On the other hand, at the moment, the peak luminosity of the radio is unlikely 
to be modified though the 'wings' of the radio spectrum are modified by the nonthermal 
electrons. Therefore, as long as we use the peak radio luminosity instead of whole structure, 
this method is not seriously defected. 
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Fig. 1. — Plot of the ratio of the radio luminosity to the 2-10 keV X-ray luminosity versus 
the X-ray luminosity. The solid lines are standard ADAF model predictions for SMBH 
masses of 10 6 — 10 9 M Q , the dotted lines are ADAF model predictions with truncations at 
25 r g , the short dashed lines and the long dashed lines are ADAF model predictions with 
winds, where p = 0.4 and p = 0.99, respectively. The open circles denote the observed 
core radio luminosity at 15 GHz, and the filled circles the spatially resolved core radio 
luminosity which are converted from ones at 5 GHz to at 15 GHz using the z/ 7//5 power law. 
The sources are spatially resolved, and the core luminosities are adopted in order to avoid 
possible contaminations due to other components such as jets. SMBH masses are denoted 
by the log scale at the top of each curve. See the text for detailed discussions. Radio flux 
values at 15 GHz of Sgr A*, NGC 1068 are quoted from Kormendy & Richstone (1995) and 
Gallimore et al. (1996), NGC 3079, NGC 3628, NGC 4151, and NGC 4388 are quoted from 
Carral et al. (1990), NGC 3031 (M81), NGC 4258, NGC 4736 (M94), and NGC 5194 (M51) 
are quoted from Turner & Ho (1994), respectively. Radio flux data at 5 GHz of NGC 1365 
and NGC 3310 are quoted from Saikia et al. (1994), and NGC 224 (M31), NGC 3377, NGC 
4374 (M84), NGC 4486 (M87), and NGC 4594 (M104) are quoted from Franceschini et al. 
(1998), respectively. X-ray fluxes (2-10 keV) of Sgr A*, NGC 224, NGC 1068, NGC 1365, 
NGC 3031, NGC 3079, NGC 3310, NGC 3628, NGC 4151, NGC 4258, NGC 4374, NGC 4388, 
NGC 4594, NGC 4736, and NGC 5194 are quoted from Baganoff et al. (2001), Trinchieri et 
al. (1999), Turner et al. (1997), Iyomoto et al. (1997), Ishisaki et al. (1996), Cagnoni et 
al. (1998), Ptak et al. (1999), Ptak et al. (1998), Nandra et al. (1997), Makishima et al. 
(1994), Colbert & Mushotzky (1999), Iwasawa et al. (1997), Forster et al. (1999), Ptak et 
al.(1998), Roberts et al.(1999), Terashima et al. (1998), respectively. 
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Fig. 2. — Core radio luminosity as a function of the SMBH mass, with SMBHs whose 
masses are available. Mass accretion rates are denoted by the log scale at each line. Lines 
correspond to same models as in Figure 1 and symbols represent same observational data 
points as in Figure 1. References for radio fluxes are shown in the caption of Figure 1. The 
mass estimates of Sgr A*, NGC 1068, NGC 3031 (M81), NGC 3079, NGC 4374 (M84), 4594 
(M104) come from Eckart & Genzel (1997), Greenhill et al. (1996), Ho (1999), Ferrarese 
& Merritt (2000), Ho (1999), Ferrarese & Merritt (2000), respectively. The estimates of 
SMBH mass in NGC 224 come from Dressier & Richstone (1988), Bacon et al. (1994). The 
mass estimate of NGC 3377 comes from Kormendy & Richstone (1995) and Magorrian et 
al. (1998). The mass of NGC 4151 have been measured by reverberation mapping method 
(Wandel et al. 1999; Gebhardt et al. 2000b). The mass estimate of NGC 4258 comes from 
Miyoshi et al. (1995). The mass estimates of NGC 4486 (M87) come from Ford et al. (1994), 
Reynolds et al. (1996), Magorrian et al. (1998), and Ho (1999). 
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Fig. 3. — Similar plot as Figure 2. But the X-ray luminosity for the 2-10 keV band is shown 
with the mass-known SMBHs. Lines correspond to same models as in Figure 1. References 
for X-ray fluxes are shown in the caption of Figure 1. References for mass estimates are 
shown in the caption of Figure 2. 



